The Dirac Quantum Cellular Automaton in one dimension: Zitterbewegung and 

scattering from potential 



Alessandro Bisioj^] Giacomo Mauro D'Ariano[j] and Alessandro Tosinj^] 
Dipartimento di Fisica delVUniversitd di Pavia, via Bassi 6, 27100 Pavia and 
Istituto Nazionale di Fisica Nucleare, Gruppo IV, via Bassi 6, 27100 Pavia 

We study the dynamical behaviour of the quantum cellular automaton of Refs. [U [2]? which 
reproduces the Dirac dynamics in the limit of small wavevectors and masses. We present analytical 
evaluations along with computer simulations, showing how the automaton exhibits typical Dirac 
dynamical features, as the Zitterbewegung and the scattering behaviour from potential that gives 
rise to the so-called Klein paradox. The motivation is to show concretely how pure processing of 
quantum information can lead to particle mechanics as an emergent feature, an issue that has been 
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I. INTRODUCTION 

The idea of reproducing the evolution of a macroscopic 
system starting from a simple rule of local interaction 
among its elementary constituents was first formalized 
in the pioneering von Neumann's paper [3 with the no- 
tion of Cellular Automaton. The automaton is a regular 
lattice of cells with a finite number of states, equipped 
with a rule that updates the cell states from time t to 
time t + 1. Such rule must be local, namely the state of 
the x cell at t + 1 depends only on the states of a finite 
number of neighboring cells at t. Cellular automata have 
been a popular topic for many years, as a new paradigm 
for complex systems, and many books have been devoted 
to the subject (see eg. Refs. [H [5]). One of the rea- 
sons of its first success, which eventually has become its 
own weakness, is the chaotic behaviour of the automaton 
dynamics [6]. 

Differently from classical cellular automata, Quan- 
tum Cellular Automata (QCA) exhibit a less chaotic be- 
haviour, which makes them predictable for large number 
of steps [7] . Here the cells are finite-dimensional quantum 
systems interacting locally and unitarily. Being local- 
ity of interactions an essential ingredient of any physical 
evolution, QCA have been considered already by Feyn- 
man as candidates for simulating physics [HI [9]. More 
recently QCA earned interest in the quantum informa- 
tion community leading to many results on its mathe- 
matical theory [10-12 , and on their general dynamical 
features [7J fl3HT6] . In quantum field theory, after the 
first appearance of a prototype of QCA in the Feynman 
chessboard [8] for solving the path-integral for the Dirac 
field, a similar framework has appeared in the work of 
Nakamura [17] motivated by a rigorous formulation of 
the Feynman path integral, and later in the seminal work 
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of Bialynicki-Birula [18] , as a lattice theory for Weyl, 
Dirac, and Maxwell fields. Then the possibility of using 
automata for describing the evolution of relativistic fields 
emerged in the context of lattice-gas simulations, espe- 
cially in the work of Meyer [19] , where a notion of "field 
automaton" first appeared, and in the papers of Yepez 
[20]. 

More recently QCA have been considered for extending 
quantum field theory [1 to the Plank scale. Similar to 
lattice-gas theories, here the quantum cell corresponds to 
the evaluation ip(x) of a quantum field on the site x of a 
lattice, with the dynamics updated in discrete time steps 
by a local unitary evolution. However, differently from 
lattice-gas theory, here the continuum limit is not taken, 
whereas, instead, the asymptotic large-scale (Fermi) evo- 
lution is considered. The main difference is then that 
Lorentz covariance holds exactly in the relativistic limit 
of small momentum, whereas generally it is distorted, in 
a fashion analogous to Refs. [2TH23] . In this context 
the one dimensional Dirac automaton has been derived 
from symmetry principles for the QCA [2] showing how 
the usual Dirac dynamics emerges at the Fermi scale, 
though relativistic covariance and other symmetries are 
violated at the Planck/ultrarelativistic scale. 

In the present paper we analyze in detail the one- 
particle sector of the automaton of Refs. [TJ [2]. Here, 
particle states are "smooth" states peaked around a mo- 
mentum eigenstate of the QCA. We will consider dynam- 
ical quantities as the particle position, momentum and 
velocity, along with their evolution both in the free case 
and in the presence of a potential, recovering typical fea- 
tures of Dirac quantum field evolution-as Zitterbewegung 
and Klein paradox-from the pure quantum information 
processing of the QCA. Recently there has been a re- 
newed interest in Dirac features in solid-state and atomic 
physics, which provide a physical hardware to simulate 
the dynamics. Zitterbewegung can be seen in the re- 
sponse of electrons to external fields [24] and can appear 
for nonrelativistic particles in a crystal [25 -27 , quasipar- 
ticles in superconductors [28] and systems with spin-orbit 
coupling [29, 30 . Proving that the oscillation behavior is 
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not unique to Dirac electrons, but rather is a generic fea- 
ture of spinor systems with linear dispersion relations, 
these works opened the way for possible simulation of 
Zitterbewegung using for example trapped ions [3TJ [32] , 
two-band crystalline structure such as graphene [33j [34] 
or semiconductors [35-39 , ultra cold atoms [40], and fi- 
nally photonic crystals [41 . On the other hand, the Klein 
paradox (tunneling of relativistic particles) provides in- 
sight in the mechanics of relativistic particles propagating 
through potential barriers, along with vacuum polariza- 
tion effects, and has been a focus in the hot topic of 
graphene as a simulator for Dirac equation, as in Ref. 
[42] , and [32] for trapped ions. Recently also microfab- 
ricated optical waveguide circuits have become an alter- 
native physical simulator for particle dynamics [43 . 
After reviewing the Dirac QCA in Id in Section [Til in 



Section [HI] we present the evolution of position and mo- 
mentum operators for the automaton, showing the Zit- 
terbewegung behaviour produced by the interference be- 
tween positive and negative frequencies. In Section IV 
we modify the QCA in order to insert a potential in the 
free evolution, and show the automaton dynamics in the 
presence of a barrier for one particle states. 

We end the paper with a summary and some conclud- 
ing remarks in SectiorfVf 



per step (c = tpjrp in dimensional units, corresponding 
to the speed of light). The quantum field can be taken 
generally as Fermionic, Bosonic, or even Anyonic. How- 
ever, in the present case it will not be relevant, since 
we will consider only single-particle states, which span 
the Hilbert space C 2 (g) ^(^), and for which we will use 
the factorized orthonormal basis where for \s) we 

consider the canonical basis corresponding to s = /,r. 
These states can be also obtained as ipl(x)\Q) upon in- 
troducing a vacuum \Q) which is annihilated by the field 
operator, and invariant under the automaton evolution. 
Similarly also Af-particle states with N > 1 can be ob- 
tained by acting with products of N evaluations of the 
field operator, building up the Fock space in the usual 
way. Notice that the evolution of the field is restricted 
to be linear, and there exists a unitary operator U such 
that the field evolution is given by ViJj s (x)V^ = Uip s (x), 
with V\Q) — whereas for product of field evalua- 
tions the evolution is given by tensor powers of U as 
V^ Sl (xi) ■ ■ .^s N (x N )V^ = U® N tp ai (xi)®. ■ .®i> SN (x N ). 

In the \s)\x) representation the unitary matrix U can 
be written as follows 



n\x — l)(x\ —im\x)(x\ 
—im\x)(x\ n\x + l)(x\J ' 



(3) 



II. THE DIRAC AUTOMATON 

The quantum automaton corresponding to the Dirac 
equation in Id, first introduced in [T], has been derived 
from the discrete automaton symmetries of parity and 
time-reversal in Ref. [2], where also the Dirac equation 
has been recovered as the large-scale relativistic limit of 
the automaton. The cell of the quantum automaton is 
given by the evaluation ip(x) of the two-component field 
operator ip, and the unitary evolution of one step of the 
automaton is given by 



ijj(x) — » Uijj(x), ip(x) 



1pr(x) 

i/ji(x) 



(1) 



where ipi and ip r denote the left and right mode of the 
field, whereas the unitary matrix U is given by 



U 



nS —im 
-im nS^ 



m 



i, 



(2) 



with S denoting the shift operator Sf(x) = f(x + 1). 
The constants n and ra in the last equation can be cho- 
sen positive. As shown in Refs. [TJ [2], the parameter 
ra plays the role of an a-dimensional inert ial mass, and 
is bounded by unit. We remark that the automaton de- 
scription is completely a-dimensional, and a conversion 
to the usual physical dimensions needs a length, a time 
and a mass, which one can take as the Planck length £p, 
the Planck time rp, and the Planck mass rap, the latter 
playing the role of the bound for the inertial mass. The 
maximal speed of propagation of information is one cell 



describing a Quantum Walk on the Hilbert space C 2 (g) 
l 2 (Z) 0. 

Tanks to the translational invar iance of U, it is conve- 
nient to move to the momentum representation 
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and U becomes 

7T 

U= I dk U{k)®\k)(kl U(k) 



k e [-7T,7T], (4) 



ne %k 

—im ne 



im 

ik 



(5) 



Notice that discreteness bound momenta to the Brillouin 
zone, as in solid-state theory. By diagonalizing the uni- 
tary matrix U(k) 



U(k)\s) k = e- is "^\s) k , 



\s)k ■= 3 



\/2 



y i - sv(k) ' 

s-s/l + sv(k) 



u(k) = arccos(n cos k) (6) 
s = ±, v(k) := dk0o(k) 



it is easy to evaluate the logarithm of U (e lH := U) as 
follows 

H = [ dk H(k)(g)\k)(k\, (7) 

J —TV 

H(k) = uj(k) (\+) k (+\ k - \-) k (-\ k ) 

= sinc _1 Co>(&)(— n sin k cr 3 + raai), 

where G{ i = 1, 2, 3 denote the usual Pauli matrices. 

The function uj(k) is the dispersion relation of the au- 
tomaton, which recovers the usual Dirac one uo(k) = 
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FIG. 1: The Dirac automaton dispersion relation in Eq. (|6| 
for four different values of the mass: m = 0.1, 0.2, 0.4, 0.8. 

\Jk 2 + m 2 in the limit k, m <C 1 and k/m^> 1 as shown 
in [2]. This is also clear in Fig. [I] where the dispersion 
relation as a function of k is reported for four different 
values of the mass. The derivative v(k) in Eq. ([6| is 
then the group velocity of the wavepacket. The s = +1 
eigenvalues correspond to positive-energy particle states, 
whereas the negative s = — 1 eigenvalues correspond to 
negative-energy anti-particle states. 

Notice that the operator H regarded as an Hamiltonian 
would interpolate the evolution to continuous time as 
U(t) = J7*, which, however, in this context should be 
considered unphysical [49] 

In the following sections we will analyze two typical 
aspects of the Dirac-field dynamics, namely the Zitter- 
bewegung and the Klein paradox. 

III. POSITION AND MOMENTUM 
OPERATORS AND ZITTERBEWEGUNG 

The QCA Q describes very precisely the Dirac field 
dynamics for customary relativistic wavevectors and en- 
ergies (consider that e.g. a ultra-high-energy cosmic ray 
has k ~ 10 -8 ) [2 . In this section we will show how effi- 
ciently it reproduces a typical feature of the one-particle 
Dirac dynamics, namely the Zitterbewegung. 

The Zitterbewegung was first recognized by 
Schrodinger in 1930 [44] who noticed that in the 
Dirac equation describing the free relativistic electron 
the velocity operator does not commute with the Dirac 
Hamiltonian: the evolution of the position operator, 
in addition to the classical motion shows a very fast 
periodic oscillation with frequency 2mc 2 and amplitude 
equal to the Compton wavelength h/mc with m the rest 
mass of the relativistic particle. This jittering motion 
first encountered in the Dirac theory of the electron was 
then shown [24] to arise from the interference of states 
corresponding to the positive and negative energies 
resulting from the Dirac equation with the trembling 
disappearing with time [45] for a wavepacked particle 
state. Zitterbewegung oscillations cannot be directly 



observed by current experimental techniques for a 
Dirac electron since the amplitude should by very small 
w 10 -12 m. However, it can be seen in a number of 
solid-state, atomic-physics, photonic-cristal and optical 
waveguide simulators, as quoted in the introduction. 

The "position" operator X providing the representa- 
tion \x) (i.e. such that = x\s)\x)) is defined as 
follows 

X = ^x(I®\x)(x\). (8) 

Generally X provides the average location of a 
wavepacket in terms of (i/j\X\i/j). The conjugated "mo- 
mentum" operator is given by 

-k{I®\k){k\). (9) 

One can verify that X and P obey the usual canon- 
ical commutation rule [X, P] = i. In the following it 
will be convenient to work with the continuous time t 
interpolating exactly the discrete automaton evolution, 
namely U l . However, all numerical results will be given 
only for discrete t, namely for repeated applications of 
the automaton unitary U in Eq. 

The time evolution of the position operator X(t) = 
U^XU 1 can be more easily computed by integrating the 
differential equation A(t) = [H, [H,X(t)]] where H was 
defined in Eq. We have 

A(t)= ( dk A(k,t)®\k)(k\ A(k,t) = e 2iH{k)t A(k) 

J —TV 

2uo 2 

A(k) = 2 — nm cos k<72 (10) 

sin a; 

which leads to 

X(t) = X(0) + Vt + Z x (t) - Z X (0) (11) 
V(k) = -v(k) 2 a 3 + vityy/l-vik) 2 *! (12) 

Z x (k,t) = -±H- 2 (k)A(k,t) (13) 

where V is the classical component of the velocity op- 
erator which, in the base diagonalizing the Hamiltonian 
([7]), is V(k) = v(k)as and is proportional to the group 
velocity v(k). Since a generic one-particle state \ip) is a 
superposition of a positive and a negative energy state, 
i.e. |^ + ) + 1^-), the evolution of the mean value of the 
position operator X(t), can be written as 

xt(t) := (mm = x+(t)+x^(t)+x™\t) 
x±(t) := (tl>±\X(0) + Vt\tl> ± ) 
xf{t) := 2K[(^+|X(0) - Z x (0) + Z x (t)ty-)] (14) 

where 5R denotes the real part. The interference be- 
tween positive and negative frequency is responsible of 
the x l ^ l {t). The magnitude of x l ^(t) is bounded by 1/m 
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FIG. 2: (Colors online) Plots of z(k) (left) and u(k)/ir (right) 
related to the oscillation amplitude and frequency of the posi- 
tion expectation value in Eq. ( 11 ). In both cases the plots are 
reported for different values of the mass (m = 0.1, 0.2, 0.4, 0.8 
from the top in the figure on the left and from the bottom in 
the figure on the right). 



(see appendix [A]) which in the usual dimensional units 
correspond to the Compton wavelength h/mc. More- 
over the stationary phase approximation shows that for 
t — »• oo the term 25ft[(^+|Zx(^)| < 0-)], which is respon- 
sible of the oscillation, goes to as 1/y/t (see ap- 
pendix |A| and only the shift contribution coming from 
2$t[(ip+]X(0) - Z x (Q)\il>-)] survives. These results show 
that x^it) is the automaton analogue of the so-called 
Zitterbewegung. As already noticed in the introduction 
this phenomenon was never observed for a free relativis- 
ts electron because of the small value of x 1 ^ t (i) which 
is bounded by the electron Compton wave length 10 -12 
m in natural units. The results of this section are in 
agreement with the one for the Hadamard walk [46] . 

In Fig. [3] we have considered the evolution of states 
with particle and antiparticle components smoothly 
peaked around some momentum eigenstate, namely 



c + |^+)+c_|^->, \^±) 



dk 



g ko (k)\±) k \k) (15) 



where c\ 



cr_ = l and g ko is a Gaussian peaked 
around the momentum ko with width a. An easy 
computation shows that for these states the shift 
contribution reduces to 25K[(^|X(0) + Z x (0)\ip)] = 
3?(<4c_)/(2tt) Jl^dk \g ko (k)\ 2 z(k) with the function 

z(k) = m cos uj(k)/ sin 2 uj(k) bounded again by the 
Compton wavelength 1/ra and the oscillation frequency 
given by uj(0)/tt (see also Fig. [2J. 



IV. EVOLUTION WITH A SQUARE 
POTENTIAL BARRIER 



In order to study the scattering with a potential, 
we modify the automaton adding a position dependent 
phase representing a square potential barrier, as in Refs. 
[46| 147] . We will provide explicitly the transmission T 
and reflection R coefficients as functions of the energy 
and mass of the incident wavepacket and of the potential 
barrier's height. We will find a general behavior inde- 
pendently on the regime, namely on the energy and mass 
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FIG. 3: Automaton evolution of a state as in Eq.( 15 ) showing 
the Zitterbewegung of the position expectation value. Top: 
m 0.15, c+ l/y/2, c_ i/y/2, k 0, and a 40" 1 . 
The calculated shift and oscillation frequency are respectively 
(ip\X(0) + Z x (0)\ip) = 3.2 and cj(0)/tt = 0.05, accordingly 
to the simulation. Middle: m = 0.15, c+ = l/>/2, c_ = 
l/\/2, &o = 0, a = 40 _1 . The calculated shift and oscillation 
frequency are and 0.13, respectively. Bottom: m — 0.13, 
c+ = y 7 ^, c_ = 1/V3, k = 10" 2 7T, a = 40" 1 . In this case 
the particle and antiparticle contribution are not balanced 
and the average position drift velocity is thus (ip+\V\ip+) + 
(ip-\V\ip-) = (|c+| 2 - \c-\ 2 )v(k ) = 0.08, corresponding to 



an average position x^(800) + (800) = 464 (see Eq. fli}). 



Notice that for t —> oo the term 23^[( / 0+|Zx(t)| / 0-), which is 
responsible of the oscillation, goes to 0. 



of the incident particle. Increasing the value of the po- 
tential barrier beyond a certain threshold a transmitted 
wave reappears and the reflection coefficient starts de- 
creasing. The width of the R = 1 region is an increasing 
function of the mass which is proportional to the gap be- 
tween positive and negative frequency eigenvalues of the 
unitary evolution. 

For a generic potential <fr(x), the unitary evolution be- 
comes 



n\x — l){x\ 
—im\x)(x\ 



—im\x)(x\ 
n\x + l)(x\ 



We will analyze the simple case <p(x) := (j)6(x) (0(x) 
is the Heaviside step function) that is a potential step 
which is for x < (region I) and has a constant value 
(f) G [0, 2tt] for x > (region II) as illustrated in Fig. [4] 
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FIG. 4: Schematic of the potential 



Let us now study the eigenvector of U<f> of the form 

\$k) = ni|+)fc|fc) + ni^|+>_ fc |fc> + TfcHn|+> Jb /|fc , > 

ni :=^2l®\x)(x\ U u :=^I®\x)(x\ 

x<0 x>0 

where Jk and k r are functions of fc. The condition 
that \<&k) is genuinely an eigenstate of U^, i.e. U^l&k) = 
e - iuj W\$ k ) : implies that 

u(k') = oj(k) - cj) (16) 
_ e~ ik y/(l + v)(l - gj - e~ ik ' y/(l - tt)(l + gj 
k ~ -e ik y/(l-v)(l-v') + e~ ik ' y/(l + v)(l + v') 
2e z ^(v cos k — i sin k) 

Ik ~ . 

-e ik y/{l- v ){l- v ')) + e ~ ik ' ^(1 + ^(1 + ^) 

with v := i>(fc) and v' := v(fc') the group velocities of the 
incident and transmitted wave. Let us now consider the 
superposition 








FIG. 5: Reflection coefficient as a function of the potential 
barrier height <j> and of the momentum k of the incident par- 
ticle state. From the top-left to the bottom-right the reflec- 
tion coefficient is depicted for different values of the mass: 
m = 0.1, 0.2, 0.4, 0.8. 



where gk (k) is a function in Cq°[— 7r,7r] which we as- 
sumed to be smoothly peaked around fco- The state at 
time t is then 

|*(t)> := J^Lg k0 (k)e-^ k ^ k ) 

and one can verify that for t <C the state is negligible 
in region II while the only appreciable contribution in 
region II comes from the term e lk ° x which describes a 
wavepacket that moves at group velocity v(ko) and hits 
the barrier form the left. When t ^> the state can 
be approximated by a superposition of a reflected and a 
transmitted wavepacket as follows 

|*(t)> ^>/3(fc ) J^ gko {k)e-^\+)_ k \k)+ 

/c\k 
— g k ,(k')e-^ k >*|+M*'> 

where we defined 
fcg s.t. uj(k f ) = uo(ko) — (/>, 

(one can check J^= Ififfc^fc')! 2 = 1)j wnose group veloc- 
ities are — v(ko) for the reflected wave packet and v(k' Q ) 
for the transmitted wave packet (see Fig. [6|. 

The probability of finding the particle in the reflected 
wavepacket is R = |/3(fco)| 2 (reflection coefficient) while 
the probability of finding the particle in the transmit- 
ted wavepacket is T = |7(fco)| 2 (trasmission coefficient). 
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FIG. 6: Group velocity of the transmitted wave packet as 
a function of the potential barrier height <j> and of the mo- 
mentum k of the incident particle state. From the top-left to 
the bottom-right the transmitted group velocity for different 
values of the mass: m — 0.1, 0.2, 0.4, 0.8. 
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FIG. 7: Reflection coefficient for m — 0.4 and momentum 
of the incident particle ko = 2 as a function of the potential 
barrier height (j) (section of plots in Figs. 5|6 for m = 0.4, 
k = 2). 



The consistency of the result can be verified by check- 
ing that R + T = 1. For k <C m <C 1 (Schroedinger 
regime) we recover the usual reflection and transmission 
coefficient for the Schroedinger equation with a potential 
step. In Fig. [5] we plot the reflection coefficient R as a 
function of (j) and k for different values of the mass m. 
Clearly when <p = we have R = and increasing <p 
while fixing k the value increases up to R = 1. One no- 
tice that when uo(k) — arccos(n) < <p < uo{k) + arccos(n) 
Eq. (16) has solution for imaginary k' which implies an 



exponential damping of the transmitted wave and pure 
reflection. By further increasing the value of <fi beyond 
the threshold uo(k) -f arccos(n), Eq. (16) have solution for 
real k' and negative uj(k'), and then a transmitted wave 
reappears and the reflection coefficient decreases. This is 
the so called "Klein paradox" which is originated by the 
presence of positive and negative frequency eigenvalues 
of the unitary evolution. The width of the R = 1 region 
is an increasing function of the mass equal to 2 arccos(n) 
which is the gap between positive and negative frequency 
solution solutions (see Fig. [I]). 

In Fig. [7] we plot the reflection R coefficient and the 
transmitted wave velocity group v(k' ) as a function of the 
potential barrier height <j) with the incident wave packet 
having ko = 2 and m — 0.4. From the figure it is clear 
that after a plateau with R = 1 the reflection coefficient 
starts decreasing for higher potentials. In Fig. [8]we show 
the scattering simulation for four increasing values of the 
potential, say (j) = 1.42, 1.55, 2, 2.4 (see the caption to 
figure for the details). 




FIG. 8: Simulations of the Dirac automaton evolution with 
a square potential barrier. Here the automaton mass is 
m — 0.2 while the barrier turns on at x = 140. In the 
simulation the incident state is a smooth state of the form 
1^(0) ) = J-^= 9k (&)|+)fc peaked around the positive energy 
eigenstate |+)fc with ko = 2 and with gk a Gaussian having 
width a = 15 _1 . The incident group velocity is v(ko) = 0.90. 
The simulation is run for four increasing values of the poten- 
tial (f). Top-Left: Potential barrier height (j) — 1.42, reflec- 
tion coefficient R — 0.25, velocity of the transmitted particle 
v(k'o) = 0.63. Top-Right: = 1.55, R = 0.75, v(k' Q ) = 0.1. 
Bottom-Left: </> = 2, R = 0.1, v(k' ) 0. Bottom-Right: 
(j) = 2.4, R = 0.50, v(k' ) = 0.33. 



tigate and visualize a wide range of fundamental pro- 
cesses. This facts, in addition to the discreteness of the 
automaton, makes of it the ideal theoretical counterpart 
for the experimental simulators in the literature. A sim- 
ilar quantum cellular automaton can be also developed 
in two dimensions [48] , corresponding to the graphene as 
quantum simulator. 



Appendix A: Bound of the oscillating term and its 
asymptotic behavior 



V. CONCLUSIONS 

In this paper we studied the dynamics of the quan- 
tum cellular automaton of Refs. [I] |2], which gives the 
Dirac dynamics as emergent in the limit of small wavevec- 
tors. We presented computer simulations and analytical 
evaluations, focusing on typical features of the Dirac dy- 
namics, in particular the Zitterbewegung and the scat- 
tering from potential. Our automaton covers all regimes 
of masses and energy-momenta, beyond the same va- 
lidity range of the Dirac equation, with the possibility 
of considering arbitrary input states, enabling to inves- 



Here we provide an upper bound for the oscillating 



term x]£ t (t) in Eq. ( 14 ) in the position operator evolution 
derived in Section ITT and we derive its behaviour for very 



long time steps. The jittering of the position expectation 
value is caused by the operator Zx(t) which in the base 
diagonalizing the automaton Hamiltonian H ^ can be 
written as 



Z x 



(t)= e 2i ^ azt Z x {k)^\k)(kl 

J — 7T 

mcosu(k) 



Z x (k) = z(k)a 2 , z(k) 



2 sin 2 uj(k) 
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with z(k) € 7r,7r) for any m/fl. By defining 
77 dfc 



|V>±> = 
we have 



/ 

J — 7T 



^±(fe)|±) fe |fe), ^(A:)GC °°[-7r,7r] 



Since, for any m ^ 0, o;(fc) has three stationary points 
in k = 0, ±tt (^(0) = ^(ivr) = and uj^(k) + 
elsewhere in the closed interval k G [— 7r,7r], with 
a/ 2 )(0), uj( 2 \±7r) 7^ 0), the stationary phase approxima- 
tion gives 



2»[(V+|Z JC (t)|V->] ^> 



/C = 0,±7T 



ig+(k)g-(k) 



7TUj( 2 \k)t 



showing that the term 23?[('0 + |Zx(t)|V ; -)], goes to as 

i/Vt. 

In order to find an upper bound for x^it) notice that 
\x^(t)\ < 2|(^+|X(0) - Z X (0) + Z x (t)\1>-)\ 

<2(|(v> + |x(o)|V-}l + \z x (o)\ + \z x (t)\) 



and, according to the expression of Zx(k), we get 
\Z x (0)\ + \Z x (t)\<2\Z x (0)\ 



\Z X (0)\< max \z(k) \=z(0) 

feG[-7T,7T] 



2m 



Now defining the Co°[— 7r,7r] test function (p(k,k f ) 
we have 



\(^\x(0M-)\ 

dS(k-k') 



v{k,k f ) 



d(Jfe - fc 7 ) 

/ u HA- 
—dfc' S(k-k')g* + (k)g.(k') 



S(k - fc') 
d 



d<p(fc,fc') 



d(fc - fc') 



d(fc - fc') 



7r(+|fc|-)fe' 



dfc 
2^ 



/(fc) := 



g*+(k)g-(k)f(k) 



/(0) 



• 2 ' 

sin 



< max |/(fc)| = /(0) 

feG[-7T,7T] 

V^l — '/Tl 2 

m 2 



which finally gives 



K\t)\<- + 



m m 



2 ' 



(Al) 
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